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We investigate three-dimensional anharmonic phonons in tetrahedral symmetry and supercon- 
ductivity mediated by these phonons. Three-dimensional anharmonic phonon spectra are calcu- 
lated directly by solving Schrodinger equation and the superconducting transition temperature is 
determined by using the theory of strong coupling superconductivity assuming an isotropic gap 
function. With increasing the third-order anharmonicity b of the tetrahedral ion potential, we find 
a crossover in the energy spectrum to a quantum tunneling regime. We obtain strongly enhanced 
transition temperatures around the crossover point. We also investigate the anharmonic effects on 
the Debye- Waller factor, the phonon spectral functions and the density profile, as a function of 
the anharmonicity b and temperature. The isomorphic first-order transition observed in KOS2O6 
is discussed in terms of the first excited state energy A, and the coupling constant A in the strong 
coupling theory of superconductivity. Our results suggest the decrease in A and increase in A below 
the first order transition temperature. We point out that the change in the oscillation amplitude 
(a- 2 ) and (xyz) characterizes this isomorphic transition. The chemical trends of the superconducting 
transition temperature, A, and A in the /3-pyrochlore compounds are also discussed. 

PACS numbers: 74.20.-z 74.25.Kc 



I. INTRODUCTION 

Recently, various low-energy properties arising from 
anharmonic ion oscillations have attracted much atten- 
tion. An ion located at the center of an over-sized 
cage oscillates with large amplitude and thus the anhar- 
monic terms in the potential energy play an important 
role. Indeed, these anharmonic oscillations are observed 
in /3-pyrochloreii~— fillcd-skuttcruditc^r— and clathratc 
compounds^— In metallic systems, such anharmonic os- 
cillations interact with conduction electrons, and due to 
their large amplitude, the electron-phonon coupling con- 
stant becomes large. The large anharmonicity and the 
strong electron-phonon coupling lead unusual relaxation 
of conduction electrons anomalous nuclear magnetic 
relaxation time ; 12 ' 13 the sound velocity anomalies jiiiirJJi 
and the strong coupling superconductivity^ 7 -"— 

These anharmonic oscillations in these systems have 
also been studied theoretically. Dahm and Ueda dis- 
cussed the anomalous temperature dependence of the 
resistivity^— and the NMR relaxation timei^ observed 
in KOS2O6 by using a single-site anharmonic phonon 
model, employing the self-consistent Gaussian approxi- 
mation for the quartic term of the ion displacement^ Re- 
cently, Yamakagc and Kuramoto extended this approach 
to the lattice problem in the same level of Gaussian 
approximation^ These works can explain the softening 
of optical-phonon frequency as temperature decreases ob- 
served in many compounds £ 1 10 i 25 ~ — 

The anharmonic ion oscillations had been discussed 
also as a possible mechanism of high transition tem- 
perature superconductivity, particularly for high T c 
cupratesi^r— As for the systems of anharmonic ion os- 
cillations in over-sized cages, importance of contributions 



of low-energy Pr ion oscillations to the superconductivity 
is pointed out in PrOs4Sbi2 £L Recently, Chang, et aZ.jSS 
discussed the superconductivity in KOS2O6 using strong 
coupling theory of s-wavc superconductivity^ also em- 
ploying the Gaussian approximation 2 ^ for anharmonicity. 

In this paper, we focus on the /3-pyrochlorc compounds 
AOS2O6 (A=K, Rb, or Cs). Monovalent A-cations are 
located at the center of the OS12O18 cages and form a 
diamond lattice structure. Lattice dynamics was investi- 
gated by neutron-scattering experiments,— ~— and a low- 
energy optical phonon is observed at about 3 meV in the 
K compound. This is basically K-cation oscillations, and 
the same phonon is also observed at around 5-7 meV in 
the Rb and Cs compounds. The root-mean-square am- 
plitude of the K oscillation turns out to be 0.12-0.14 A at 
zero temperature from the elastic neutron scatterin g 25 ' 26 
while much smaller values are reported for Rb and Cs 
compounds^ 5 - 

The first-principle calculations indicate that the K ion 
potential has large anharmonicity and is very shallow 
along [111] and three other equivalent directions^ It is 
also calculated that the first excited state has the excita- 
tion energy of ~8 K and K-cation oscillation amplitude 
is as large as 1 A at zero temperature. These values 
are quite different from the experimental data, indicat- 
ing that some parameters in the ion potential are not so 
realistic. In the present paper, we will first systematically 
analyze the effects of anharmonicity of the ion potential 
on the ion dynamics in KOs206- 

The /3-pyrochlore compounds reveal superconducting 
phase transition and the transition temperature T c 's are 
9.6 K, 6.3 K and 3.3 K for A=K, Rb, and Cs, respec- 
tively. These values of T c are inversely related to the 
A-cation size, i.e., T c is the highest for the smallest ion: 
potassium. The symmetry of the gap function is consid- 
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ered to be fully gapped s-wave i 12 i 18 ' 19 It has been also 
shown that the electron-phonon coupling is largei 3 ' 18 ' 20 
Thus, it is expected that the conduction electrons on the 
cages strongly interact with the anharmonic oscillations 
of ion inside the cage in these systems. 

For KOs 2 Og, in addition to superconducting transi- 
tion, there exists a first-order structural transition at 
T p = 7.5 K. At this transition, no sign of symme- 
try breaking has been observed^ The oscillation of lo- 
cation seems less anharmonic below T p as indicated by 
electrical resistivity, specific heat jump at T c in magnetic 
fields and the mean free path estimated from the upper 
critical field iJ c2 .— Recently, we have proposed that this 
is driven by a sudden change in K-cation oscillation am- 
plitude driven by intersite ion interactions^ The ampli- 
tude of the low-energy excited state with xyz symmetry 
jumps at T p , but this does not change the Tj, point-group 
symmetry. 

The main purpose of this paper is to clarify the proper- 
ties of anharmonic oscillations in tetrahcdral symmetry, 
which is the local symmetry for ^4-cations in AOs 2 Og, 
and the superconductivity mediated by these anharmonic 
oscillations. In order to fully take into account the 
anharmonicity and anisotropy, we will solve the three- 
dimensional Schrodinger equation for an anharmonic os- 
cillator in the tetrahcdral symmetry. Using these exact 
phonon cigenstates, wc will then discuss the strong cou- 
pling superconductivity assuming an s-wave pairing. 

This paper is organized as follows. In Sec. UH we will 
calculate the energy spectrum of the anharmonic poten- 
tial problem and discuss various thermodynamic and dy- 
namical quantities and their dependence on temperature. 
Anharmonic effects on Debye- Waller factor will also be 
discussed. Section [III] is devoted to the discussions for 
strong coupling superconductivity mediated by the an- 
harmonic ion oscillations discussed in Sec. |TTJ In Sec. 
IIVI wc will discuss the relevance of the present results to 
the /3-pyrochlorc compounds. We also apply our theory 
to discuss the changes in phonon dynamics at the first- 
order transition at T p . Finally Sec. [V]is a summary of 
this paper. 
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where R = (X, Y, Z) is the real-space displacement of 
the ion from the equilibrium position. M is the mass of 
the ion. Throughout this paper we set M/m e = 71748, 
where m e is the mass of electron, except for discussions 
in Sec. IIV A[ and this corresponds to the mass of K 
ion. Oo and B are coefficients of second- and third-order 
terms in the ion potential, while C\ and C 2 are isotropic 
and cubic fourth-order terms, and wc ignore the higher- 
order potential terms of 0(R 5 ). To study dependence on 
potential parameters, it is useful to rewrite Hamiltonian 
(fTJ) into a dimensionless form by renormalizing displace- 
ment and energy by their units. As for the energy unit, 
we choose the energy of harmonic phonon corresponding 
the second-order term htto. The unit of the length is 
chosen as ao = %/MVLq = aB\/2(e 2 /2a.B)(m e /Af)/M2o- 
Here, as is the Bohr radius as = h 2 /m e e 2 ~ 0.53 A. 
In these units, Hamiltonian (fTJ) is transformed to 
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(3) 



V(r) = -w |rf + bxyz + ci|r| 4 + c 2 (x 4 + y 4 + z 4 )(4) 

where r = (x, y, z) = R/ao. b and c\ and c 2 are di- 
mensionless constants and ujq is introduced to vary the 
second-order term for later purpose, but ljq = 1 for the 
main part of this paper. In Sees. IIIII and ITVl we will use 
different values of ujq to examine the effects of changes in 
the second-order term^ Throughout this paper we set 
%Q,q = 44 K which corresponds to the energy scale of 
the optical-phonon frequency observed in /3-pyrochlore 
compounds ^ 25 ' 27 and thus ao ~ Ob/vIO. 

In the limit of b = c\ = c 2 = 0, Hamiltonian ([3]) can be 
diagonalizcd by using creation and annihilation operators 
as, 



II. ANHARMONIC PHONONS 
A. Model 

In this paper, we investigate an anharmonic oscilla- 
tion of K ion in KOs 2 Oe. Our model is anharmonic local 
phonons at each lattice point and we assume the local 
symmetry is tetrahedral one which corresponds to the 
case of K-site symmetry in KOs 2 06- In tetrahedral sym- 
metry, in addition to spherical and cubic fourth order 
terms, there exists a third-order anharmonic term which 
breaks inversion symmetry and the Hamiltonian is given 
by 



H = a\a x + (Aa y + a\a z + 
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where 



= 7f (^m + - t m) witn A* = 



(5) 



y or z. 



Eigenstatcs are labeled by three occupation numbers as 
H\n x ,n y ,n z ) = (n x + n v + n z + 3/2)\n x ,n y ,n z ) and 
a^a^nxjUy,^} = n^ln^, n y , n z ). 

For nonzero &, c%, and c 2 , we diagonalizc Hamiltonian 
(U) numerically in the restricted Hilbert space spanned 
by {\n x , n y ,n z )} with n x + n y + n z < n max . In this paper, 
we use n max = 40 which corresponds to the Hilbert space 
with 12341 states and check the convergence by compar- 
ing the results for n max = 50 including 23426 states. We 
employ this approach rather than the conventional self- 
consistent Gaussian approximatio n 23 ! 24 . This is because 
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the conventional approximation does not work for the po- 
tential, Eq. (UJ since the third-order term cannot be de- 
coupled as in the fourth-order terms. It is also important 
that, as will be shown later, energy differences between 
adjacent eigenstate multiplets are not the same, and this 
property cannot be described by the self-consistent Gaus- 
sian approximation. 

B. Energy spectrum 

In Fig. [TJ we show a low-energy part of the energy 
spectra of Hamiltonian @ as a function of the third- 
order anharmonicity b for c\ = 0.04 and c-i = 0.01. The 
ground state is always nodeless in the space r and there- 
fore singlet, while the first excited states are triplet, cor- 
responding to s and p orbitals for the case of isotropic 
potential. As b increases, the energy of the lowest singlet 
excited state (hereafter, this state will be referred to as 
s' state) decreases and shows an anticrossing with the 
ground state at b = b* ~ 1.7. This kind of anti-crossing 
behavior does not occur in the one-dimensional anhar- 
monic potential problem: V(X) = aX 2 + bX 4 l&22. In- 
terestingly, near b = b* the five lowest-energy states are 
well separated from other states in the energy spectra. 
Thus, the validity of our previous five-state toy model is 
justified around this crossover region^ 

For b > b* , the excited triplet states, which are p-wave 
like, are nearly degenerate with the ground state. This 
means that these four states are localized away from the 
origin in the four directions: [111], [HI], [111], and [111]. 
For the illustration, we show the potential form along 
[100], [110], and [111] directions in FigH It is clear that 
an off-center potential local minimum emerges along [111] 
direction as b increases. The energy of the first excited 
state from the ground state A is small, owing to smallness 
of the quantum tunneling probability between different 
valleys. These four states form sp 3 orbitals and their en- 
ergy can be well described by considering the quantum 
tunneling of the ion between the four potential minimum 
positions. Hereafter, we call these states in this parame- 
ter region (b > b*) as quantum tunneling states. 

C. Observables 

We now investigate the temperature dependence of the 
oscillation amplitude of these anharmonic ions. Ther- 
modynamic average of observable A is calculated by the 
formula: 

(A) = J2™n(n\A\n), (6) 

n 

where \n) is the eigen state of Hamiltonian ((TJ) with 
eigenenergy E n and w„ is its Boltzmann weight w n = 
exp(—E n /T)/Y^ m ex.p(—E m /T), and T is the tempera- 
ture. In the summation in Eq. (JB|), we discard states 
with E n > 0.8n max wo, since the high-energy parts of our 
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FIG. 1: (Color online) Energy eigenvalues vs the third-order 
anharmonic parameter b for ci = 0.04 and C2 = 0.01. s, p, 
d, and / denote the approximate orbital symmetry, which is 
exact in the isotropic harmonic potential case. 
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FIG. 2: (Color online) Potential along [100], [110], and [111] 
directions for various b. a = 0.04 and C2 = 0.01. r = rn, 
where n being the unit vector along each direction. 



energy spectrum is not correct due to the cutoff n max 
and the high-energy part does not matter as long as the 
low-temperature properties of the system arc concerned. 

In Fig |3l we show the temperature dependence of fluc- 
tuations (a; 2 ) = (r 2 )/3 and (xyz). Note that since xyz 
is invariant quantity in tetrahedral symmetry, this does 
not vanish except for 6 = 0. As the temperature de- 
creases, both (x 2 ) and \(xyz)\ decrease and saturate to 
finite values for b < b*. On the other hand, for b > b* , 
these quantities increase at low temperatures suggesting 
the quantum tunneling state. 

To see the change with b more clearly, we calculate the 
ion density p(r) = J2 n w n \^> n (r)\ 2 at position r. Here 
^„(r) is the wave function at position r which is normal- 
ized as / d 3 x\^ n (r)\ 2 = 1. In Figs. IUa)-|3|d), we show 
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FIG. 3: (Color online) Temperature dependence of thermo- 
dynamic average (a) (a; 2 ) and (b) (xyz). ci = 0.04 and 
c 2 = 0.01. 



p(r) along [111] direction, r = r(l, 1, l)/y/3, for several 
temperatures. p(r) in other directions has monotonic 
temperature dependence and is peaked at r = (0, 0, 0) as 
in Fig. Ufa). As b increases, a part of the weight of p(r) 
shifts to the position around r ~ —4. Interestingly near 
b = b* [Fig. life)], the temperature dependence of p(r) 
at r r~j —4, is non- monotonic, i.e., as T decreases, p{— 4) 
increases first, but decreases below T ~ 0.3 ~ A. In Fig. 
Efd), the p(0) is suppressed at low temperatures reflect- 
ing the nature of the quantum tunneling states, since the 
potential minimum away from the origin is deep enough 
as shown in Fig. [2] 
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FIG. 4: (Color online) Atomic density profile p(r) along 
r(l, 1, l)/\/3 for various b. ci = 0.04 and c 2 = 0.01. 

D. Debye- Waller factor 



Now we discuss the Debye- Waller factor, which con- 
tains information about the amplitude of ion oscilla- 
tion and it is observed by elastic neutron-scattering 
experiments^ Consider an atom whose position is rep- 
resented by d in a unit cell and denote its displacement 
by Ud = (uj, u d , u z A ). Then, the Debye- Waller factor at 
the scattering wave vector Q for this atom e~ Wdi - c ^ is 



given by 

cxp[-VF d (Q)] = (expHQ • u d )). 
In the harmonic approximation, Eq. ([7]) reduces to 

^f(Q) = ^((Q-u d ) 2 ). 
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Using W / d(Q), the static structure factor F(Q) is repre- 
sented as 



F(Q) = ^6 d exp(iQ-d)cxp[-Wd(Q)] 



(9) 



where & d is the averaged scattering length of the atom at 
position d. 

When anharmonicity is not negligible, Eq. (jHJ is not 
sufficient and we must take into account higher-order 
terms. Up to the third order in (Q • Ud), we obtain 



wf } (Q) = §((Q-u d )^ 



l 
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((Q • u d ) 3 



(10) 



As an example, let us consider K atoms in KOS2O6, which 
constitute a diamond sub-lattice. There are two sites in 
the unit cell, d^ = (0,0,0) and d B = § (1,1,1). Here a 
is the lattice constant. From the symmetry arguments, 
we obtain 



^((Q-u d ) 2 ) = i|Q| 2 (|u d | 2 )^lU d , 
i((Q-u d ) 3 ) = Q x Q y Q z {u x A u y A u A ) = 



(11) 
(12) 



In the diamond lattice structure, W' A (Q) = W' A (Q) = 
W'(Q) and W d ' A (Q) = -W£ B (Q) = W"(Q). Using Eqs. 
(ITT|) and (fT2"|) and setting b = b dA = & df3 , we obtain 

|F(Q)| 2 = 26 2 exp[-2W"(Q)] 



1 + cos 



Q • (d A - d B ) - 2W"(Q) h3) 



Thus, the effective tU(Q) becomes 



Weff(Q) = W"(Q)--log2 



- log 



COS 



-Q • (d A - dfl) - W"(Q) 



(14) 



It is important to note that since the W"(Q) depends on 
the wave vector as Q x QyQz, the effect of anharmonicity 
is anisotropic. For example, at Qm = (1, 1, l)27r/a, W e s 
is given by 



Eq. (fTU), |^(Q)| 2 = at Q = (2,0,0)27r/a and also 
at (2, 2, 2)2n/a due to the interference factor 1 + cos[Q • 
(dA — de)] when W"(Q) is set to be zero. For b 7^ 0, 
(u A u A u A ) is finite and this leads to nonvanishing W"(Q) 
at Q = (2,2,2)27r/a but W"(Q) = at Q = (200)27r/a, 
since W"(Q) is proportional to QxQyQz- Therefore, 
|i^(Q)| 2 is nonvanishing at Q = (2,2,2)27r/a while re- 
mains zero at Q = (2, 0, 0)27r/a. 

We show the calculated W e ff(Qin) for several values 
of b as a function of temperature in Fig. [5] In Fig. [5ja) , 
the first term in Eq. (fT5|) is shown and W e ff(Qin) f° r 
b > and b < are shown in (b) and (c), respectively. It 
is easily seen that, for small b, W c ff(Qin) in Figs. EJb) 
and [5{c) are essentially the same as W'(Qm) in (a), 
while, for larger b, Wofr(Qni) is enhanced (suppressed) 
for b > in (b) [b < in (c)]. For other directions 
of Q with Q x QyQ z = 0, W e ff(Q) does not include the 
contribution of the third order term W"(Q) and thus, 

lUeff = W. 

Figure [6] shows the wave- vector dependence of |F(Q)| 2 
for several temperatures and b. Note that F(Q) at 
Q = (2,0,0)27r/a vanishes due to the extinction rule 
mentioned above, while the intensity at Q = (2,2,2)27r/a 
appears for nonzero b although its strength is weak. Con- 
cerning 6 dependence, ^(Q)! 2 for b = 1.74 is much 
suppressed for large |Q| and low temperatures, which 
is originated from the large displacement for b > b* . 
As for the temperature dependence, |F(Q)| 2 increases 
as T decreases for b — and 1.58, while decreases for 
b = 1.74 > b*. This temperature dependence of |-F(Q)| 2 
for b > b* is also understood from the temperature de- 
pendence of W e ff in Fig. El where W e g for larger b > b* 
increases as T decreases. 

We note that by careful analysis of the neutron data, 
as we have demonstrated the effects of the third-order 
term on the Debye- Waller factor, we can obtain useful 
information about anharmonicity in the real materials. 
In KOS2O6, since there are a lot of kinds of atoms such 
as the oxygen and osmium, our results presented in this 
section cannot directly be compared with the neutron 
data of KOS206- It is desired to carry out high reso- 
lution neutron scattering experiments using the single 
crystal and analyze the results by taking into account 
the anharmonicity. 



Woff(Qlll) = Weff(-Qlll), 



W(Q m ) - - log il + sin [2W" (Qm) 



whereas at Q200 = (2, 0, 0)27r/a, we obtain 

4. 



Wc ff (Q 



200; 



W'(Q 
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E. Ion dynamics 

Now let us investigate the dynamics of the anharmonic 
' ion oscillation. In Fig. [3 we show the phonon density of 
(15) states (DOS) F(u), which is calculated from the correla- 
tion function D, 



(16) 



-tW'M = y^(w„ - w m ) 



(n|a;, x |m)(m|x Ai /|n) 



ui — E, n 



E„ 



rij 



(17) 



The third-order contribution (u A u A u A ) also modifies The DOS is given by its imaginary part, F(uj) = 
the extinction rule of the structure factor. According to — ImD xx (u))/Tr. Notice D xx (ui) = D yy {ui) = D zz {ui) in 
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FIG. 5: (Color online) Temperature dependence of the expo- 
nent of the Debye-Waller factor for ci = 0.04 and C2 = 0.01. 
(a) f'(Q m ). (b) Weff(Qiu) for b > 0. (c) VK eff (Qm) for 
b < 0. Inset in (c): schematic potential form along the [111] 
direction in (b) and (c) for large b cases. 



the tetrahedral symmetry. We set the phenomenological 
broadening parameter rj = 0.26. Although one can try 
more sophisticated analyses by introducing dissipation 
by electron-phonon couplings or couplings with other de- 
grees of freedom to calculate the imaginary part of the 
self-energy ri(uj) f 3r ^ 8 it is sufficient at this stage to re- 
strict ourselves in the phenomenological level as long as 
characteristic properties of the DOS are concerned. 

Figure [3a) shows the phonon DOS F(lu) for the case 
b = and the eigenenergies relative to the ground state 
energy indicated by small bars. It is clearly seen that the 
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FIG. 6: (Color online) Q dependence of |F(Q)| 2 for ci = 0.04 
and C2 = 0.01 at temperatures (a) T — 4.0, (b) T = 1.0, and 
(c) T = 0.02. For each of Q, |F(Q)| 2 values for 6 = 0, 1.58, 
and 1.74 are shown from left to right. For Q = 2n(2, 0, 0)/a, 
F(Q)| 2 is always zero. 



position of the peak decreases as the temperature de- 
creases and this softening is consistent with the study 
based on self-consistent Gaussian approximations > 23 i 24 
The single peak in F{uj) is actually constituted of many 
Lorentzian peaks each of which has a width of r\ = 0.26. 
In this case of b = 0, the transition matrix elements are 
sizable only between neighboring "multiplets" of the cor- 
responding harmonic system. Because of the fourth-order 
anharmonicity, c\ and C2, the degeneracy of each mul- 
tiplet is only approximate and slightly lifted, and also 
the energy separation between multiplets increases with 
increasing energy. Since the distribution of Lorentzian 
peak positions is quite smooth and the energy depen- 
dence of Boltzmann weight is monotonic, they finally 
form a single broad peak. Its position shifts from a low 
energy at low T where only the ground state is an im- 
portant initial state to a higher energy at higher T where 
higher excited states contribute more importantly. 

For b > 0, qualitatively different behaviors appear as 
shown in Figs. EJb)-[7Jf). The transition matrix elements 
are now finite for many more pairs of eigenstates (n|a;|m) 
for nonzero b. This is because the point-group symmetry 
changes from cubic Oh to tetrahedral Td and the eigen- 
states |n) do not have definite parity related to inversion 
symmetry r — > — r. In this point-group symmetry, for 
example, x and yz are bases of the same Ti irreducible 
representation. As a result, the matrix element such as 
(yz\x\s) is nonvanishing (here "s" means an isotropic s- 
wave state, i.e., A\ in Td point group). In particular, the 
ground state has transition matrix elements much larger 
than in the case of b = 0. Thus, even at the lowest tem- 
perature, more than one peaks are present in F(u)) for 
b > 0, while the 6 = case has no additional visible 
high-energy peak in F(lu) as shown in Fig. [TJa), since 
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FIG. 7: (Color online) Phonon DOS F(ui) for various fe's fixing 
ci = 0.04 and c 2 = 0.01. The data correspond to T =10, 9.2, 
8.4, 7.6, 6.8, 6.0, 5.2, 4.4, 3.6, 2.8, 2.0, 1.8, 1.6, 1.4, 1.2, 1.0, 
0.8, 0.6, 0.4, 0.2, and 0.1, from bottom to top and each curve is 
shifted up by 0.005. Energy eigenvalues relative to the energy 
of the ground state are shown by short bars in each panel. 



the ground state has nonvanishing matrix elements only 
with the odd-parity excited states and their magnitude 
is small. 

For clarity, let us concentrate on a few specific peaks 
to understand their origin for 6^0. F(uj) for the lowest 
temperature T = 0.1 in Fig. [3c) shows four peaks below 
lu ~ 3. As it is easily checked by comparing the energy 
eigenvalues indicated by bars in Fig. [7J the positions of 
these peaks correspond to the excitation energies. Ac- 
tually, the lowest three peaks in Fig. [7J (c) correspond 
to the transitions to the Ti states which can be traced 
back to the p-wave sates in the one-phonon multiplct, 
the d-wave states in the two-phonon multiplet, and the 
p- or /-wave states in the three-phonon multiplet in the 
harmonic oscillator (see Fig. [T] where the label of the 
multiplets in the harmonic oscillator are indicated). 

There are also peaks or shoulder-like structures origi- 
nating from the transition between the s' and T2 states 
for b > 0. Among them, the transition between the first 
excited p-wave states to s' state is clearly seen in Figs. 
[Jjb) and[7]Jc). Since the energy difference between these 
excited states A' is smaller than that between the ground 
state and the first excited state A, these shoulder- like 
structures appear at energies lower than the main peak 
position which corresponds to the transition from the 
ground state and the first excited states for the inter- 



mediate temperatures. In Fig. [Z^d), A' is so small that 
the structure at w ~ A' is not clearly visible but the 
slope F(uj)/uj near ui ~ is enhanced. 

In Figs. EJe) and[7jf), the first excited states are al- 
most degenerate with the ground state, so that the low- 
est energy peak in the figures does not correspond to the 
transition to the first excited states but that between the 
first excited states and s' state. 



III. STRONG COUPLING 
SUPERCONDUCTIVITY 

In this section, we will discuss the strong coupling 
theory of superconductivit y 39 ! 40 in which the attractive 
force between electrons is mediated by the anharmonic 
phonons discussed in the previous section. The main 
issue of this section is the effect of the phonon anhar- 
monicity on the superconducting transition temperature 
T c . 

We assume that the phonon DOS is given by F(oj) 
obtained in Sec. Ill El i.e., we ignore the renormaliza- 
tion due to the electron-phonon coupling. Generally, 
electron-phonon couplings lead to frequency rcnormaliza- 
tion and a finite life time of phonons, both of which are 
described in the renormalization of F(u>). This renormal- 
ization is interesting and generally should be included, 
but the present calculation reproduces correct qualitative 
behaviors and we leave full self-consistent calculation as 
a future problem. It is noted that the frequency renor- 
malization is taken into account such that the potential 
parameters are chosen to reproduce a renormalized fre- 
quency. The superconducting transition temperature T c 
is obtained by applying the conventional strong coupling 
theory of s-wave pairin g 39 ! 40 to our anharmonic phonon 
system coupled to isotropic electron gas. 

Throughout this section, we use Matsubara formalism, 
which is efficient to determine T c . Matsubara formulation 
has also advantage that the phcnomenological broaden- 
ing factor 77 used in Sec. Ill El is not necessary. 



A. Gap equations 

Following the conventional theory of strong coupling 
superconductivity— the transition temperature T c for 
isotropic s-wave gap Age is determined by the gap equa- 
tion, 

A S c(ie m ) = -ct 2 T^ K(ie n )D(ie m - ie„)A sc (ie n ),(18) 

n 

/oo 
d£G(£,ie n )G(£,-ie n ), (19) 
-OO 

where the electron Green's function is given by 

G «-*-> = fr.-fUfc.) - < 2o) 
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and the normal self-energy is given by 

/oo 
d£G(Z,ie n - iv m )D{iv m ). (21) 

Here D(iv m ) = D^^{w m ) is the phonon Green's func- 
tion. e n = (2n + l)nT and v m = 2rmrT are the fermionic 
and bosonic Matsubara frequencies, respectively, a 2 is 
proportional to the square of the electron-phonon cou- 
pling constant times electron DOS at Fermi energy. In 
our model, it is a frequency-independent quantity and 
is set as a 2 = 595 K/A 2 , which corresponds to 0.38 in 
our units of energy and length. As in the conventional 
theor y 39 i 40 , the normal self-energy, Eq. (|2ip is essen- 
tially given by the second-order perturbation theory and 
K(ie n ) is analytically given by^ 2 - 
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To determine T c , we define 



FIG. 8: (Color online) Superconducting transition temper- 
ature T c vs the phonon energy A for the harmonic case 
(22) b = ci = C2 = 0. The dashed line corresponds to T C M given 
by Eq. 



27|) and the solid line represents T c given by Eq. 



M mn = —a 2 T^K(ie m )D(ie m — ie n ) y/ K (ie n ) , (23) and the characteristic energy scale uj\ og is given by 



4' n = A sc {ie n )\/ K ( ie n), 
and solve eigenvalue problem numerically 

A(T)?/' m = J2 Mmn ^ n ~ 



(24) 
(25) 



(26) 



2a 2 f 00 F(w) 
logwiog = — — / duj logo;, 



(29) 



2a 2 v 



A ^ — ' _EV r) — -E/U 

E m >E„ 



|(n|x|m)| 2 log (E m -E n )(30) 



The Allen-Dynes formula is given by 



With decreasing temperature, A(T) increases and T c is 
given by the condition: A(T C ) = 1. 



B. Superconducting transition temperature in 
harmonic systems 

Before discussing the effects of phonon anharmonicity 
on T c , we show in Fig. [8]T C in the harmonic system as 
a function of the first excited state energy A = luq. The 
dash line is T c calculated by McMillan's formula^, which 
is valid in the weak-coupling regime, while the solid line is 
that by Allen and Dynes,— which gives a good estimation 
of T c in the extremely strong coupling regime. Let us 
summarize two analytic formulae for the two opposite 
limits. The McMillan's formula is given by 



T, 



AD 



o.isvA^ 2 



(31) 



T C M = 0.8w log exp ( - 



1 + A 
A 



(27) 



where the Coulomb pseudo-potential is not taken into 
account. Here, the dimensionless coupling constant A is 
given by 

\ = a doj =a — |(n|x|m)| , (28) 



where (a; 2 ) = 2a 2 / °° ujF(uj)duj/ A. 

In the case of harmonic phonons, A can be analytically 
calculated and is given by A = a 2 /A 2 .— This means the 
smaller A is, the larger A is realized. The McMillan's for- 
mula (|2"T|) indeed reproduces T c for large A as expected, 
while for small A, T c approaches the value of the Allen- 
Dynes formula (011) (= 0.18V0.38 ~ 0.11) as shown in 
Fig. [5] Note that for the harmonic phonons, (u 2 ) = A 2 . 
This relation is valid only for harmonic oscillators as will 
be discussed below. 



C. Superconducting transition temperature in 
anharmonic systems 

Now, we discuss the effects of anharmonicity of phonon 
dynamics on T c . Figure [9] shows the dependence of T c on 
the third-order anharmonic term b of the ion potential, 
Eq. Q for c 2 = 0, 0.01, and 0.03, at fixed a = 0.04. 
All the three cases exhibit a pronounced peak around 
the crossover value b = b* as discussed in Sec. |TT] for 
each C2, and T c 's are strongly suppressed in the quan- 
tum tunneling states, i.e., for large b > b* . The solid 
lines in the figure represent the McMillan's formula T C M 
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FIG. 9: (Color online) T c vs b for ci = 0.04 and c 2 =0, 0.01, 
and 0.03. The solid lines show T c 's of McMillan's formula 
(|27[1 . 6*'s for each of c 2 are indicated by arrows. Inset: zoom 
up in the large b part. 

[Eq. (f2T|)]. The agreement with the McMillan's formula 
indicates that T c is qualitatively given by the conven- 
tional theory of strong coupling superconductivity. Note 
that A and LU\ og vary with the temperature in the pres- 
ence of anharmonicity and T C M are calculated by using 
these values evaluated at T = T c . For b > b* , however, 
the discrepancy becomes large as shown in the inset of 
Fig. [9j and the effects of anharmonicity and the nature 
of the extreme strong coupling regime appear, which will 
be discussed later in Sec. IIII El 
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FIG. 10: (Color online) A, wi og) and A vs 6 for c 2 = 0, 0.01, 
and 0.03. The open symbols represent A and the hlled sym- 
bols with a solid line represent ui\ os . Lines without symbols 
represent A. A and u;i og are calculated at T — T c . 
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D. Competition between the energy scale and 
coupling constant 

In Sec. IIII Dl we discuss the origin of the peak struc- 
ture based on the McMillan formula^ 

The McMillan's formula is expressed with the two pa- 
rameters, A and wi og , as shown in Eq. (|27[) . In Fig. 
I10[ we show the b dependence of A and a>i og calculated 
at T c together with A. It is noticeable that A shows a 
steep increase above b* and uji og is almost the same as A 
for b < b* . From these facts, one can easily understand 
that T c calculated by Eq. (|zT7| decreases in the quan- 
tum tunneling states. The peak structure in T c in Fig. 
[9] is a result of the competition between the suppression 
of wi g(A) and the enhancement of A. Thus, the peak 
structure is realized at b ~ b* , where A is large enough 
and simultaneously wi og is not vanishingly small. 



E. Strong coupling limit 

In an extreme strong coupling regime A > 1, Allen 
and Dynes showed that T c is approximately given by Eq. 
([51]) rather than the McMillan's formula ||27)l44 When 



FIG. 11: (Color online) A, A' vs A for c 2 = 0.01 and five sets 
of ljo and ci. The open (filled) symbols represent A' (A). A 
is controlled by varying b. 



Eq. (|3ip is applied to our local phonon problem, we 
need care about the quantity X(uj 2 ). As pointed out by 
Hardy and Flocken^ this quantity is related to the f- 
sum rule and this implies T c is universal in any potential 
and depends only on the mass of the ion and a 2 as shown 
in Fig. [5J However, this is not true in our anharmonic 
phonon system. From numerical calculations, instead of 
(a; 2 ), we find that T c for b 3> b* is approximately given 
by 

T c = 0.18A\/A. (32) 

We note that the origin of large A for b > b* is the strong 
suppression of A for b > b* and other high- energy states 
do not play an important role. Thus, yj (ui 2 ) overestimate 
the energy scale of T c . These facts s uppo rt the use of A 
instead of the averaged frequency \J (co 2 ). 
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FIG. 12: (Color online) A and A s ; vs A for c- 2 = 0.01 and five 
sets of wo and Ci. The open (filled) symbols represent A (A s /). 



F. Additional channel of attractive interactions 

As we noted before, there is a crossover at b = b* to 
the quantum tunneling states and this means there are 
five states at low-energy region of the spectra as shown 
in Fig. [T] In the harmonic case, A is basically deter- 
mined by the matrix element (n|x|0) between the ground 
state |0) and the threefold-degenerate first excited states, 
and the corresponding excitation energy A. Since four 
low-energy excited states are nearly degenerate with the 
ground state around b = b* , the singlet excited state \s') 
also has a noticeable matrix element (n|x|s'). This pro- 
vides another channel of attractive interaction and con- 
tributes to T c . It is important to note that this additional 
contribution never appears in the one-dimensional model 
and is a direct consequence of the crossover of the ground 
state at b = b* . 

In Fig. [TTJ wc show A and the energy difference be- 
tween s' state and the p-wave like first excited states: 
A' = E s i — E p . They are calculated with varying b, and 
plotted as a function of the energy difference A between 
the first excited state and the ground state. A increases 
roughly as ~ A~ 3 / 4 for small A. A' shows minimum 
around A ^ 0.2 and this reflects the anti-crossing in the 
energy spectra. Corresponding to this, the coupling con- 
stant \ s > related to s' state shows a maximum around 
A ~ 0.1 as shown in Fig. [12J Here, we define \ s > as 

A s , = 2a 2 J2(w„ - w s ,)\(n\x\s')\ 2 /(E s , - E n ), (33) 

n—p 

where the n summation is taken over the three first ex- 
cited states. Note that the smaller A' is, the larger A s / 
is obtained. Depending on the magnitude of A S ' , total A 
has a small bump around the A ~ 0.1 as a function of A. 
This observation clearly shows that there is a new chan- 
nel of interaction via the excitations between the first 
excited states and the s' state. 



Deviations from Allen-Dynes formula near 

b = b* 



As we discussed in Sec. IIII F[ the enhancement of A s ' 
influences the total A and also affects T c itself. In Sec. 
IIII G[ we derive a formula which includes the contribution 
of A s ' . In order to take into account this contribution in 
Eq. (f3"2"j). we assume the DOS given by 



F(w) 



X p S(uj - A) + X S 'S(uj - A') 



(34) 



Since the enhancement of A s < takes place in a finite but 
extremely strong coupling regime, the total A is large and 
we follow the discussion by Allen and Dynes to estimate 
the lower bound of T c , by assuming a trial gap function 
Asc(^„) = A§ c <$| Wre | iirT i£ Using the form ([3i]). we ob- 
tain the lower bound of T r as 



87T 2 T 2 



A P A 2 + (A.. - 1)A' 2 



[A P A 2 + (A., - 1)A' 2 ] 2 + 4(A P + A S ,)A 2 A' 



(35) 



In the limit of A s ' — > Eq. (|35j) reduces to Eq. 
T c = A^/Ap"/(27r). Although the more sophisticated gap 
function can give the almost correct coefficient 0.18 in- 
stead of (27r) -1 ,— it is enough to consider the simplest 
trial gap function in order to examine the contributions 
of the s' state to T c . In Eq. (|35[) . we assume X p (^- 1) 
corresponds to the coupling constant originating from 
the transition between the ground state and p-wave like 
states (first excited states) and \ s > <C A p is the part from 
s' state [Eq. (|33l) ]. Since A p ^> X s >, we can obtain a 
simpler approximate form of Eq. (|35[) . This is obtained 
by simply replacing A in Eq. (|32p by A e ff, 



A ff — A p 



A s , 



l + 2Ap £ 



(36) 



Figure [T31 shows the A dependence of T c /A for two typ- 
ical sets of parameters in log-log scale. As shown with 
the dashed line, the Allen-Dynes formula shows a straight 
line. For small A, the calculated values agree with the 
McMillan's formula, while for extremely large A they ap- 
proach the line of the Allen-Dynes formula. For the inter- 
mediate A > 1, however, the results show a bending and 
this bending is larger in Fig. fT3Tb) than in Fig. HUT a). 
Our two-channel formula (|35[) can describe the bending 
as indicated by the solid line. Although the lower bound 
of T Ci Eq. ([33]) is not in very good agreement with the 
calculated T c , it captures the overall behavior of T c for 
A > 1. It is also noted that, although the McMillan's 
formula (f27|) shows an increase in T c /A for large A, the 
increase is too steep. The present two-component anal- 
ysis confirms the existence of the additional channel of 
the interaction near the crossover point in this system. 
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FIG. 13: (Color online) Various T c formulas fl27J, (Ell), (EU, 
and (|36p vs A. Formula (|27[1 with o;i og replaced by A is also 
plotted for comparison, (a) cjq — 0.6, ci = 0.04, and C2 
O.f) (b) cj = 1, ci = 0.036, and c 2 = 0.01. 

IV. DISCUSSIONS 

Let us now discuss the implications of our calculations 
to compare with characteristic properties experimentally 
observed in /3-pyrochlore compounds. One point is the 
chemical trends among the three member compounds, 
and we discuss why the potassium compound has the 
highest T c of superconductivity. Another point is about 
the question why the superconducting phase is not so 
much affected by the isomorphic structure transition at 
T p in the phase diagram. We will also propose a possible 
change in the K-oscillation profile at T p and the effects 
of the transition on superconductivity. 

A. Chemical trends in /3-pyrochlore compounds 

Among the three /3-pyrochlore compounds AOS2O6 
(A=K, Rb or Cs), the K compound has the highest 
T C K =9.6 K of superconductivity and the strongest an- 
harmonicity in the A-cation oscillation dynamics as ob- 
served in the neutron-scattering experiments: the Dcbye- 
Wallcr factor of K ion is much smaller than of Rb and Cs 
compounds 2 ' 5 and the softening of the low-energy phonon 
peak is also the strongest in KC^Og^ 7 - The Rb com- 
pound has the next strongest anharmonicity and the sec- 
ond highest T c Rb =6.3 K, while the Cs compound has the 
weakest anharmonicity and the lowest T c Cs =3.3 K. The 
ratio of T c is approximately 3:2:1. Thus, the anharmonic- 



ity in the ion dynamics and the value of T c are related. 
Let us examine this point in our results and also check 
if one can explain, at least qualitatively, the trends of 
other important quantities, electron mass enhancement 
and phonon energy. 

A crucial difference among the three compounds is the 
size of A-cation; the K ion has the smallest size, followed 
by Rb, and Cs is the biggest ion. Since the size of the 
surrounding OS12O18 cage essentially does not change, 
the K ion has the largest space inside the cage,— leading 
to strongly anharmonic oscillations. This reflects in dif- 
ferent shapes of the A-cation potential, as shown in the 
calculation by Kunes, et al.^- although the point-group 
symmetry is common. 

In our theory, the compound-dependent potential 
shapes are modeled by adjusting parameters in the po- 
tential, Eq. @ or (01 . In the following, we will show that 
T c values and the ion dynamics in the three compounds 
are naturally explained by appropriate choice of poten- 
tial parameters. Throughout Sec. IIV A[ we will use some 
parameters explicitly shown with physical dimensions if 
necessary. 

Let us first consider the trend of the value b/b* among 
the three compounds, recalling the anharmonicity ap- 
pearing in the Debye- Waller factor and the phonon en- 
ergy observed in the neutron experiments. 25 ' 27 As shown 
in Fig. [31 the anharmonicity grows with b/b* . This im- 
plies that the b/b* value is the smallest for Cs, then for 
Rb, and the largest for K. Note that the largest ratio 
is still smaller than or at most equal to the crossover 
value b/b* = 1, since the observed K-cation density pro- 
file does not change its peak position from the equilib- 
rium position ! 46 ' 47 

Secondly, we consider the trend of the second-order po- 
tential parameter uq. Figure [Ml shows the b dependence 
of T c for several values of luq. The first-principle calcu- 
lation of the ion potential clearly shows that ujq is very 
small for K while larger for Rb and the largest for Cs4£ 
This trend is consistent with the fact that the K com- 
pound has the strongest anharmonicity and the highest 
T c as far as b < b* , as shown in Fig. [TJ] It is also impor- 
tant to note that the peak position ~ b* becomes smaller 
as wo decreases. 

By keeping these variations in b/b* and ujq in mind, 
let us discuss our results of A, A at T = T c and T c itself 
for each of the three compounds. In order to make dis- 
cussion simple, we assume the electron phonon coupling 
and the electron density of states are the same for the 
three compounds and a 2 = 1897 K/A 2 . The atomic mass 
of each ^4-cation is Mk = 71748m e , = 156839m e , 
and M Cs = 243890m e , respectively. We have adjusted 
potential parameters, Eq. for each compound to 

reproduce the phonon energy and the effective ion os- 
cillation variance (x 2 ) cxp determined from the Debye- 
Waller factor 2^22. As shown in Sec. IIID| Wd(Q) gener- 
ally contains the contribution of the third-order fluctua- 
tions (xyz) in addition to (x 2 ). However, since the ex- 
perimental data is the Q-averaged value of Wd(Q)/|Q| 2 , 
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FIG. 14: (Color online) T c vs b for five values of ujo fixing 
ci = 0.04 and c 2 = 0.01. 



it is sufficient to consider the second-order fluctuations 
(W d (Q)/\Q\ 2 ) = i(|u d | 2 ) = ^(a; 2 )- Since our potential, 
Eq. or dU, includes four parameters, we fix B = 9324 
K/A 3 and Ci = 4C 2 = 3332 K/A 4 for simplicity and vary 
It is noted that f^o is expected to become smaller as 
the size of the alkali cation decreases, and, for fixed B, 
b/b* becomes smaller as fio increases. 

Tabic |T] shows the list of the basic quantities, A, A, and 
T c for the three compounds. Let us first discuss T c . It is 
quite sure that one needs also to include a high-energy 
phonon, which plays a role to increase the T c by about 3- 
5 K for all the three compounds, to obtain T c consistent 

(2) 

with the experimental results. T c is one calculated with 
including an additional phonon with energy hfl^ = 260 
K. We have set the corresponding dimensionless coupling 
constant A/j = 0.256 in the calculations. The results of 

(2) 

Tc arc quantitatively consistent with the experimental 
values T C K = 9.6 K, T c Rb = 6.3 K, and T c Cs = 3.3 K. Here, 
we can reproduce the experimental values of T c by assum- 
ing the same a 2 , Xh and flh for all the three compounds 
and we conclude that the origin of the difference in T c 
for the three /3-pyrochlore compounds is due to the dif- 
ference mainly in the anharmonicity of the alkali cation 
oscillations^ 

Now, let us discuss the electron effective mass. In our 
theory, the electron mass enhancement factor is related to 
Atot as 1/ztot = 1 + A to t- A's obtained from experiments 
and previous theoretical studies are A K =1.6-2. 4 ) 3 i 17 i 20 i 22 
for K, A Rb =1.0-1.3^32 f or Rb and A Cs = 0.78,^ for Cs, 
respectively. Our estimation of A to t = A + A/, = A + 0.256 
with A in Table U is qualitatively consistent with these 
values, although we have not fit T c for each compound. 
The experimental values of specific-heat coefficient 7 is 
7k = 70, 7Rb = 45, and 7c s =41 mJ/mol K 2 while the 
band calculations predicted 70 ~ 10 mJ/mol K 2 for all 
the three compounds ! 20 ! 50 Thus, the experimental mass 



TABLE I: Local potential parameters and (a; 2 ) at T = 0, 
A at T = T c , A, T c , and T c (2) for the three /3-pyrochlore 
compounds. The other parameters are B = 9324 K/A , 
C 2 = Ci/4 = 3332 K/A 4 , and a 2 = 1897 K/A 2 . 

hQ (K) b/b* (x 2 )(r = 0)(A 2 ) A A(K) T C (K) T C (2) (K) 

~K 26T 058 0.0152 1.47 38.5 6~45 WlT~ 
Rb 54.6 0.28 0.0050 0.34 56.2 0.82 5.74 
Cs 74.8 0.11 0.0024 0.12 75.2 <0.03 3.37 



enhancement factor 7k/7o in KOS2O6 is 1.5-1.7 times 
larger than in RbOs206 and CsOs20g 5 i In our calcu- 
lation, zS/z t K ot ~ 1.71 and z^J zf ot ~ 1.98, which are 
semi-quantitatively consistent with the experimental val- 
ues. From our results, the enhanced mass enhancement 
in KOS2O6 is attributed to the proximity to the crossover 
to the quantum tunneling state and also the large oscilla- 
tion amplitude, which effectively enhances the electron- 
phonon coupling. 

In order to obtain the correct value of the mass en- 
hancement (~ 7 for KOS2O6 and ~ 4 for the other two), 
we should also take into account the electron-electron in- 
teractions, but this is beyond the scope of the present 
study and we leave it for a future problem. Nevertheless, 
it is important to note that A is significantly enhanced 
near b = b* by anharmonic oscillations and that this is 
the central reason why the mass enhancement in KOS2O6 
is much larger than in Rb and Cs compounds. 

As for the phonon frequency A, optical modes related 
to the A-cation oscillations are observed at energy around 
5-7 meV at room temperature in the inelastic neutron 
scatterin g 25 ' 27 and their energies show strong softening 
as the temperature decreases. This corresponds to our 
results shown in Fig. [7] For KOS2O6, the phonon en- 
ergy at T = 1.5 K is about 3.1 mcV, which is quite small 
compared with the value of 5.5 meV at T = 300 Kj2! Fur- 
thermore, the low-energy phonon peaks in KOS2O6 are 
rather broad compared with Rb and Cs compounds^ 7 - 
which suggests the proximity to the crossover point. The 
specific-heat experiment 3 - and the photo-emission spec- 
troscopy data 5 ^ also support these results. 

In our calculation for KOS2O6 with parameters in Ta- 
ble [H the peak position in the phonon spectrum is about 
70 K ~ 6 meV at T = 300 K and this shifts to A = 38.5 
K as T decreases. This means the potential parameters 
in Table U can reproduce not only T c and A but also the 
temperature dependence of the spectra. 

As for the chemical trends, we have calculated su- 
perconducting transition temperature for the three /?- 
pyrochlore compounds adjusting parameters in the po- 
tential, Eq. @ , and we can reproduce the phonon energy 
and amplitude related to the average Debye- Waller fac- 
tor. We have obtained quantitatively consistent values of 
T c with the experimental ones, using the same electron- 
phonon coupling, the electron DOS and also the same ad- 
ditional high-energy phonon for all the three compounds. 
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FIG. 15: Schematic temperature-magnetic field phase dia- 
gram taken from Ref. 3. The solid and dotted lines represent 
H C 2 and T v , respectively. The dashed line indicates the ex- 
trapolation of low-T H C 2 to the higher-temperature region. 




FIG. 16: (Color online) ujq dependence of characteristic quan- 
tities, when T c is fixed to 0.05 by tuning b. (a) A, A, and b. 
(b) The second and third-order moments of ion oscillation 
(x 2 ), and \(xyz}\. ci = 0.04 and c 2 = 0.01. 



B. Changes at T v 



As mentioned in Sec. HI KOS2O6 exhibits an isomor- 
phic first-order transition at T p = 7.5 K.— i 4 *^ 7 - In Sec. 
IIV B| let us investigate the effects of this transition on 
the phonon dynamics, using a set of constraints offered 
by the experimental results. We also discuss its effects 
on superconductivity. 

Electric resistivity shows a concave temperature 
dependence at high temperatures, and this is at- 
tributed to the strong coupling to phonons with strong 
anharmonicity^^ At the magnetic field H = 14 T, the 
resistivity is suddenly suppressed by 25 % at T p , and 
shows a different T-dependence proportional to T 2 at 
T c < T < T p .— This indicates that the electron-phonon 
scattering processes are reduced significantly at X^2. This 
is also consistent with the reduction in the specific-heat 
jump at T c in magnetic fields above 8 T.— 

In the magnetic field-temperature phase diagram, T p 
remains essentially insensitive to magnetic field H as 
shown in Fig. 1151 The upper critical magnetic field H& 
is suppressed below T p , but extrapolating to the region 
above T p , it seems that it also vanishes at the position 
very close to T C (H = 0). 

This fact cannot be explained if the characteristic en- 
ergy scale, wiog or A is common across T p : the extrap- 
olated H C 2 line should vanish below T C (H = 0) since A 
is reduced below T p as discussed above. Therefore, the 
above fact implies the enhancement of the characteristic 
energy scale A or wi og below T p . Indeed, Chang, et al.^ 
assumed a slightly increased Einstein energy to fit the 
specific heat data and our previous study also predicted 
the increase in the oscillation energy^ 

Let us discuss the changes in A, A, | (x 2 ) | , and | (xyz) | at 
T p based on our results, modeling the isomorphic transi- 
tion by a sudden change in the two potential parameters, 



u>o and b. This change corresponds to the variation in 
the mean field part of inter-site ion interactions, which 
was discussed in Ref. 34 and also changes in the oxygen 
positions and the lattice constant 4£ Two parameters are 
chosen under the constraints (i) two different parameter 
sets above and below T p lead to the same T c , and (ii) 
A is smaller below T p , both of which are the experimen- 
tal constraints, and we also assume (iii) no significant 
change in the electron band structure across the tran- 
sition and thus a 2 is also unchanged, and (iv) b < b* 
for both below and above T p implied by the result of 
the electron-density profile obtained by the x-ray and the 
neutron experiment s 46 ' 47 as discussed before. In the fol- 
lowing, we discuss the change across T p based on these 
constraints. 

In Fig. 1161 characteristic quantities are shown for fixed 
ci and C2 with varying ujq and simultaneously b such that 
those give the same T c = 0.05. Since b < b* which is the 
assumption (iv) , T c is well approximated by the McMillan 
formula (|27|) . in which T c is determined by the two factors 
wiog ~ A and A. 

Firstly, it is noted that A increases as loq increases. 
One might expect a suppression of A as ujq increases, but 
the key point is that we simultaneously tune b to fix T c 
unchanged and this means larger b is necessary at larger 
wq. This increase in b overcomes the competing effect 
of the increase in ojq, and finally A increases. Secondly, 
the energy of the first excited state A increases as wo 
decreases. This is naturally understood by noting that 
A becomes large for smaller b. Thus, from the constraint 
(ii) the above results imply that the phonon energy is 
enhanced below T p and, indeed, this is consistent with 
the previous theoretical studies ! 22 ' 34 

Most interestingly, the second- and third-order fluctu- 
ations (x 2 ) and I (xyz) I change oppositely as shown in 
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Fig. [TBTbh As ujq decreases, (a; 2 ) increases for the most 
of the luq range, while \(xyz)\ decreases monotonically. 
This is because small wo corresponds to small b as ex- 
plained above. These results indicate that (a; 2 ) is slightly 
enhanced below T p while | (xyz) | is suppressed. 

In order to illustrate these anisotropic fluctuations, the 
density isosurface for p(r) = 0.05 and the density (color) 
map are shown for two different parameters at T — 0.002 
in Fig. [T7J for (a) and (b), uj a = 0.8 and 1.0, and b = 
1.44 and 1.66, respectively, and these two sets give the 
same T c = 0.045, and we choose the two parameter sets 
to emphasize the change in the oscillation profile. As 
expected from the values of b, the case (a) has the smaller 
A = 0.647 and the larger A = 0.662, while A = 0.782 and 
A = 0.429 for the case (b). 

One can see that the density isosurface shows more 
anisotropic character in the case (b), where the "spikes" 
sticking out along four [111] directions are sharper than 
those in the case (a). This aspect reflects in the val- 
ues of (x 2 ) and (xyz). Indeed, (x 2 ) at T — T c is larger 
(x 2 ) = 0.645 in (a) than 0.626 in (b), while \(xyz)\ at 
T = T c is smaller \(xyz)\ — 1.93 in (a) than 2.35 in 
(b). Figures EO^-EIe) show p(r) along [111], [001], and 
[110] directions, respectively, for the same two data. It is 
clearly seen that |(xj/.z)| in (b) is larger than that in (a) 
as shown in (c), and (a; 2 ) in (a) is larger than in (b) as 
shown in (d) and (e). 

These changes in the anisotropy of density distribu- 
tion is a characteristic nature of the first-order isomor- 
phic transition at T p . This transition is isomorphic, con- 
sistent with the experiments^ because the amplitude of 
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FIG. 17: (Color online) Ion density isosurface for p(r) = 0.05 
and the density map on the section with (110) plane at T = 
0.002. (a) Lu = 0.8, b = 1.44 and (b) lu = 1, b = 1.66. (c) 
p(r) along [111] direction corresponding to (a) and (b), (d) 
along [001] and (e) along [110] direction. p(— |r|)'s in (d) and 
(e) are not shown since p(r) = p(— r) along these directions. 



anisotropy changes but it does not break the point group, 
translational, and any other symmetries. Our results 
suggest that (x 2 ) is slightly enhanced while \(xyz)\ is 
suppressed below T p in KOS2O6. Recent high-resolution 
neutron-scattering experiment shows the increase in (a; 2 ) 
across T p as T decreases^ which is consistent with our 
result. 

Furthermore, although it is naively expected that an 
increase in (a; 2 ) corresponds to the enhancement of A, 
this does not necessarily hold in this system, since the 
value of A is sensitive to the value of 6, especially near 
the crossover point b* . Our results show that the value 
of I (xyz) I plays more important role for the enhancement 
of A than (a: 2 ). In order to detect this change in the den- 
sity distribution across T p , it is important to perform de- 
tailed neutron-scattering experiments and to analyze the 
results with taking into account the third-order term in 
Debye- Waller factor as discussed in Sec. Ill Dl It is worth- 
while to examine the Q dependence of the Debye- Waller 
factor to extract the anisotropy of the oscillations, since 
Q-averaged Debye- Waller factor hinders the anisotropy 
and anharmonicity. 

Finally, we comment on the possibility of re-entrant su- 
perconducting transition. Experimentally, Hiroi et a/., 3 
observed a re-entrant superconducting transition around 
H = 7 T as shown in Fig. [151 This can be understood 
from H C 2(T = 0) in the strong coupling theory of super- 
conductivity: 



H c2 (T = 0) oc T 2 (l + A) 2 (l + 1.44 



(37) 



for the clean limits In our analysis, the A decreases and 
wiog increases below T p , while the T c is approximately the 
same for the potential parameters above and below T p . 
Thus, H C 2(T = 0) for the low-T parameter set is smaller 
than that for the high-T set. Therefore, we expect a 
H C 2(T) curve similar to that shown in Fig. [15] This 
result does not change when we use the formula for the 
dirty limit^ 



V. SUMMARY 

We have investigated in the present paper anharmonic 
phonons and strong coupling superconductivity in f3- 
pyrochlorc compounds, AOS2O6 (A=K, Rb, or Cs). 

First, we have solved the Schrddingcr equation of 
the three-dimensional anharmonic phonon in tetrahedral 
symmetry. The main issue is the importance of the third- 
order anharmonicity bxyz in the ion potential allowed for 
this symmetry. We have determined the energy spectrum 
of the anharmonic phonon as a function of the third- 
order anharmonicity b. We have found that there exists 
a crossover of the ground state to the quantum tunneling 
state for b > b* . We have pointed out non-monotonic 
temperature dependence of ion density profile near the 
crossover point b*. 
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Secondly, we have calculated the transition tempera- 
ture of superconductivity T c mediated by these anhar- 
monic phonons. We have found that the enhancement 
of T c near 6 = 6*. Its 6 dependence can be well fitted 
by the McMillan formula at 6 < 6* . For 6 > 6* we have 
found that there exists an additional channel of pairing 
interaction, which turns out to originate from low-energy 
excited states appearing at 6 ~ 6* . We have analyzed its 
contribution and derived an approximated formula of T c 
in the strong coupling limit. 

We have also discussed the chemical trends of T c in 
the /3-pyrochlore family AOs 2 O e (A=K, Rb, or Cs). The 
main difference among the three members are different 
values of 6/6* and the second-order term in the ion po- 
tential wo- By assuming the same electron-phonon cou- 
pling constant a 2 and the same high-energy phonon for 
all the three compounds, the differences in T c and the 
energy of the first excited phonon states A have been 
quantitatively explained only by the difference in the lo- 
cal anharmonic potential. 

Finally, we have discussed the effect of the first-order 
transition observed in KC^Og- The changes in the den- 
sity distribution of K cation at the first-order transition 
have been discussed based on the experimental data ob- 
tained so far. Especially, we have found that (x 2 ) and 
change differently across T p . Our results sug- 
gest the increase in the first excited phonon energy A 



and the reduction in the dimcnsionless electron-phonon 
coupling constant A across T p . While the latter is con- 
sistent with the experimental results, the former has not 
been observed and this is our prediction for the experi- 
ments. In order to detect the increase in A, it is impor- 
tant to carry out the high-resolution inelastic neutron 
scattering experiments. We hope our results shed lights 
in /3-pyrochlore compounds both on the strong coupling 
superconductivity and on the anharmonic oscillations ob- 
served. 



Acknowledgments 

The authors would thank T. Dahm, K. Ueda, J. Ya- 
maura and Z. Hiroi for grateful discussions. They also 
acknowledge the international workshop "New Develop- 
ments in Theory of Superconductivity" held at Insti- 
tute for Solid State Physics, University of Tokyo, dur- 
ing June 22- July 10, 2009 for providing an opportunity 
of discussion with other participants. This work is sup- 
ported by KAKENHI (Grants No. 19052003 and No. 
20740189) and also by the Next Generation Super Com- 
puting Project, Nanoscience Program, from the MEXT 
of Japan. 



* Electronic address: hattori@issp.u-tokyo.ac.jp 

1 Z. Hiroi, S. Yonezawa, and Y. Muraoka, J. Phys. Soc. Jpn. 
73, 1651 (2004). 

2 Z. Hiroi, S. Yonezawa, J. Yamaura, T. Muramatsu and Y. 
Muraoka, J. Phys. Soc. Jpn. 74, 1682 (2005). 

3 Z. Hiroi, S. Yonezawa, Y. Nagao and J. Yamaura, Phys. 
Rev. B 76, 014523 (2007). 

4 T. Hasegawa, Y. Takasu, N. Ogita, M. Udagawa, J. I. Ya- 
maura, Y. Nagao, and Z. Hiroi, Phys. Rev. B 77, 064303 
(2008). 

5 T. Goto, Y. Nemoto, K. Sakai, T. Yamaguchi, M. Akatsu, 
T. Yanagisawa, H. Hazama, K. Onuki, H. Sugawara, and 
H. Sato, Phys. Rev. B 69, 180511(R) (2004). 

6 K. Kaneko, N. Metoki, H. Kimura, Y. Noda, T. D. Mat- 
suda, and M. Kohgi, J. Phys. Soc. Jpn. 78, 074710 (2009). 

7 H. Kotegawa, H. Hidaka, T. C. Kobayashi, D. Kikuchi, 
H. Sugawara, and H. Sato, Phys. Rev. Lett. 99, 156408 
(2007). 

8 B. C. Sales, B. C. Chakoumakos, R. Jin, J. R. Thompson, 
and D. Mandrus, Phys. Rev. B 63, 245113 (2001). 

9 I. Zerec, V. Keppens, M. A. McGuire, D. Mandrus, B. 
C. Sales, and P. Thalmeier, Phys. Rev. Lett. 92, 185502 
(2004). 

H. Matsumoto, T. Mori, K. Iwamoto, S. Goshima, S. 
Kushibiki, and N. Toyota, Phys. Rev. B 79, 214306 (2009). 

1 Y. Kasahara, Y. Shimono, T. Shibauchi, Y. Matsuda, S. 
Yonezawa, Y. Muraoka, and Z. Hiroi, Phys. Rev. Lett. 96, 
247004 (2006). 

2 M. Yoshida, K. Arai, R. Kaido, M. Takigawa, S. Yonezawa, 
Y. Muraoka, and Z. Hiroi, Phys. Rev. Lett. 98, 197002 



(2007). 

Y. Nakai, K. Ishida, H. Sugawara, D. Kikuchi, and H. Sato, 
Phys. Rev. B 77, 041 101 (R) (2008). 

Y. Nemoto, T. Yamaguchi, T. Horino, M. Akatsu, T. 

Yanagisawa, T. Goto, O. Suzuki, A. Donni, and T. Ko- 

matsubara, Phys. Rev. B 68, 184109 (2003). 

I. Ishii, T. Fujita, I. Mori, H. Sugawara, M. Yoshizawa, K. 

Takegahara, and T. Suzuki, J. Phys. Soc. Jpn. 78, 084601 

(2009). 

Y. Nakanishi and M. Yoshizawa (private communications). 

M. Briihwiler, S. M. Kazakov, J. Karpinski, and B. Bat- 

logg, Phys. Rev. B 73, 094518 (2006). 

I. Bonalde, R. Ribeiro, W. Bramer-Escamilla, J. Yamaura, 

Y. Nagao, and Z. Hiroi, Phys. Rev. Lett. 98, 227003 

(2007). 

Y. Shimono, T. Shibauchi, Y. Kasahara, T. Kato, K. 

Hashimoto, Y. Matsuda, J. Yamaura, Y. Nagao, and Z. 

Hiroi, Phys. Rev. Lett. 98, 257004 (2007). 

Y. Nagao, J. Yamaura, H. Ogusu, Y. Okamoto, and Z. 

Hiroi, J. Phys. Soc. Jpn. 78, 064702 (2009). 

J. Chang, I. Eremin, P. Thalmeier, and P. Fulde, Phys. 

Rev. B 76, 220510(R) (2007). 

J. Chang, I. Elemin, and P. Thalmeier, New J. Phys. 11, 
055068 (2009). 

T. Dahm and K. Ueda, Phys. Rev. Lett. 99, 187003 (2007). 
A. Yamakage and Y. Kuramoto, J. Phys. Soc. Jpn. 78, 
064602 (2009). 

K. Sasai, K. Hirota, Y. Nagao, S. Yonezawa, and Z. Hiroi, 
J. Phys. Soc. Jpn. 76, 104603 (2007). 

R. Galati, C. Simon, P. F. Henry, and M. T. Weller, Phys. 



16 



Rev. B 77, 104523 (2008). 

27 H. Mutka, M. M. Koza, M. R. Johnson, Z. Hiroi, J. I. 
Yamaura, and Y. Nagao, Phys. Rev. B 78, 104307 (2008). 

28 N. M. Plakida, V. L. Aksenov, and S. L. Drechsler, Euro- 
phys. Lett. 4, 1309 (1987). 

29 J. R. Hardy and J. W. Flocken, Phys. Rev. Lett. 60, 2191 
(1988). 

30 V. H. Crespi and M. L. Cohen, Phys. Rev. B. 48, 398 
(1993). 

31 G. D. Mahan and J. O. Sofo, Phys. Rev. B 47, 8050 (1993). 

32 S. Manalo, H. Michor, G. Hilscher, M. Briihwiler and B. 
Batlogg, Phys. Rev. B 73, 224520 (2006). 

33 J. Kunes, T. Jeong, and W. E. Pickett, Phys. Rev. B 70, 
174510 (2004); J. Kunes and W. E. Pickett, Physica B 
378-380, 898 (2006). 

34 K. Hattori and H. Tsunetsugu, J. Phys. Soc. Jpn 78, 
013603 (2009). 

35 In those cases of ujo 7^ 1, we keep the same energy unit of 
44 K, instead of renormalization of energy. This simplifies 
our discussions about too dependence in later sections. 

36 S. W. Lovesey, Theory of Neutron Scattering from Con- 
densed Matter (Clarendon Press, Oxford, 1984), Vol. 1.. 

37 K. Hattori and K. Miyake, J. Phys. Soc. Jpn. 76, 094603 
(2007). 

38 M. Takechi and K. Ueda, J. Phys. Soc. Jpn. 78, 024604 
(2009). 

39 G. M. Eliashberg, Zh. Eksp. Teor. Fiz. 38, 966 (1960) [Sov. 
Phys. JETP 11, 696 (I960)]. 

40 D. J. Scalapino, J. R. Schrieffer, and J. W. Wilkins, Phys. 
Rev. 148, 263 (1966). 

41 A. A. Abrikosov, L. P. Gorkov, and I. Ye. Dzyaloshinskii, 
Quantum Field Theoretical Methods in Satistical Physics, 
2nd ed. (Pergamon Press, Oxford, 1965). 

42 G. Bergmann and D. Rainer, Z. Phys. 263, 59 (1973). 

43 W. L. McMillan, Phys. Rev. 167, 331 (1968). 

44 P. B. Allen and R. C. Dynes, Phys. Rev. B 12, 905 (1975). 



45 When A is expressed by the quantities with the dimension, 
A = (2a 2 /A)(fi/2MA), which is, indeed, dimensionless. 

46 J. Yamaura, S. Yonezawa, Y. Muraoka, and Z. Hiroi, J. 
Solid State Chem. 179, 336 (2006). 

47 K. Sasai, M. Kofu, R. M. Ibberson, K. Hirota, J. Yamaura, 
Z. Hiroi, and O. Yamamuro, J. Phys.: Condens. Matter 22, 
015403 (2010). 

48 The result in Ref. 33 suggests a slightly negative value of 
u)o for KOS2O6. However, the first excited state then has an 
energy too small compared with the experimental estimate 
from the specific heat data (Ref. 3) and the neutron scat- 
tering (Refs. 25 and 27). Correspondingly, calculated {x 2 ) 
is too large compared with the experimental results and 
its temperature dependence is qualitatively different from 
the experimental one (Refs. 25 and 26). In this paper, we 
consider the cases of o>o > for simplicity. We hope that 
more detailed information on the ion potential is provided 
by the neutron scattering experiments in the future. 

49 In this point whether the additional phonon is acoustic 
or optical mode does not make an essential difference. To 
improve estimation of T c , we need the detailed information 
of each compound on the electron-phonon coupling a 2 , Qh, 
Xh , and other parameters. 

50 Z. Hiroi, J. Yamaura, S. Yonezawa, and H. Harima, Phys- 
ica C 460-462, 20 (2007). 

Note that this value of the mass enhancement factor in- 
cludes not only the contributions of electron-phonon inter- 
actions but also those of electron-electron interactions. 

52 T. Shimojima, Y. Shibata, K. Ishizaka, T. Kiss, T. Tagoshi, 
A. Chainani, T. Yokoya, X.-Y. Wang, C. T. Chen, S. 
Watanabe, J. Yamaura, S. Yonezawa, Y. Muraoka, Z. Hi- 
roi, T. Saitoh and S. Shin, Phys. Rev. Lett. 99, 117003 
(2007). 

53 J. P. Carbotte, Rev. Mod. Phys. 62, 1027 (1990). 



